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ABSTRACT 

The extreme outer Galaxy (EOG), the region with a Galactic radius of more 
than 18kpc, is known to have very low metallicity, about one-tenth that of the 
solar neighborhood. We obtained deep near-infrared (NIR) images of two very 
young (~0.5Myr) star-forming clusters that are one of the most distant embed- 
ded clusters in the EOG. We find that in both clusters the fraction of stars with 
NIR excess, which originates from the circumstellar dust disk at radii of <0.1 AU, 
is significantly lower than those in the solar neighborhood. Our results suggest 
that most stars forming in the low-metallicity environment experience disk dis- 
persal at an earlier stage (<lMyr) than those forming in the solar metallicity 
environment (as much as ~5-6Myr). Such rapid disk dispersal may make the 
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formation of planets difficult, and the shorter disk lifetime with lower metallicity 
could contribute to the strong metallicity dependence of the well-known "planet- 
metallicity correlation", which states the probability of a star hosting a planet 
increases steeply with stellar metallicity. The reason for the rapid disk dispersal 
could be increase of the mass accretion rate and/or the effective far-ultraviolet 
photoevaporation due to the low extinction; however, another unknown mecha- 
nism for the EOG environment could be contributing significantly. 

Subject headings: infrared: stars - planetary systems: protoplanetary disks 
- stars: pre-main-sequence - open clusters and associations: individual (Digel 
Cloud 2) - stars: formation 



1. INTRODUCTION 



Becaus e the EOG is known t o have an extreme environment such as very low metallicity, 
[0/H]~— 1 ([Rudolph et al.ll2006l ). it can be studied as a "laboratory" f or the star formation 



process in an environm ent totally different from the solar neighborhood ([Ferguson et al. 



2006, 



Kobayashi et al.ll2008l ). As a part of our studies of star formation in the EOG (jYasui et al. 



1998 



20081 ). we explore in this paper the lifetime of circumstellar disks in a low- metallicity 



environment. 



The frequency of disks harboring stars within a young embedded cluster, the disk frac- 
tion, has been studied for many clusters of various ages as an indicator of the disk lifetime 
(IHaisch et al.ll2001al : iHernandez et al.l 120071 ). The most common means of estimating the 
disk fraction is the detection of infrared emission from warm and hot dust (typically a 



few 100 to 1500 K) in 



(Kenvon fc Hartmann 



1999 



;he in ner part of the disk (within a stellocentric distance of ~1AU) 



19871 ) with wavele ngths of 2-24 /xm u sing A-band (2.2 /xm) (ILada 



Hillenbrand! I2005T). L-band (3.5 /xm) ( Haisch et al.ll2001al) . and Spitzer IRAC 3.6-8 /xm 



and/or MIPS 24 /xm (ISicilia-Aguilar et al.ll2006l ) filters. The disk fraction was found to de- 
crease as a function of age, showing that the inner disk lifetime is about 5-10 Myr for nearby 
embedded clusters with the solar metallicity. Because (sub-)mm continuum observations of 
cold dust ( ~10K) in the outer disk (> 50AU) correlate very well with observations of inner 
dust disks ([Andrews fc Williams! 120051 ). the entire disk, from 0.1 to 100 AU, is thought to be 



1 Based on data collected at Subaru Telescope, which is operated by the National Astronomical Observa- 
tory of Japan. 
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dispersed almost simultaneously (within 
is also thought to be 5-10 Myrs. 



40 5 yr). Therefore, the lifetime of the entire disk 



The disk lifetime in a low-metallicity environment is of special interest because of the 
discovery that the probability of a star hosting; a planet increases wi t h stellar metallicit y 
("planet-metallicity correlation") (I Gonzalez! 1 19971 ; iFischer et al.l 120051 ; ISantos et al.ll2003l ). 
However, the disk lifetime in an environment with lower metallicity than the solar neighbor- 
hood has not been studied. By directly studying the disk lifetime in the EOG, it may be 
possible to obtain clues to the mechanism of this "planet-metallicity co rrelation." Although 



the disk fraction has been studied for a low-metallicity ([0/H]~— 0.5) (lArnault et al 



198 



star-forming cluster, 30 Do radus, in the Large Magellanic Cloud (LMC) (IMaercker fc Burton 
20051 ; IMaercker et al.ll2006l ). the target stars were limited to only very massive ones (~20 M Q ) 
because of the limited sensitivity and spatial resolution for this extragalactic cluster. Be- 
cause the EOG is much closer than nearby dwarf galaxies, it allows us to obtain accurate disk 
fraction s in a low-metallicity en vironment to ~0.3 M Q , which is below the characteristic peak 
of IMF (jEfmegreen et al.ll2008l ). This mass detection limit corresponds to the l owest mass of 



the t arget stars for the planet-metallicity correlation (e.g., FGK-type stars) (iFischer et al. 
20051 ) and is therefore critical for understanding the correlation. 



2. OBSERVATION 



Our targets, the Digel Cloud 2 clusters (jYasui et al.l 120081 ; iKobayashi et al.l 120081 1 are 
among the most distant embedded clusters in the EOG at R g ~ 19kpc. They are located 
i n a giant molecu lar cloud, Cloud 2, one of the most distant molecular clouds in the EOG 
(jPigel et al.lll994l ). The meta llicity [O/H] of Cloud 2 is measured at — 0.7dex from the ra- 
dio molecular emission lines fjLubowich et al.l 120041 ) . which is comparable to th at of LMC 
(~-0.5 dex) and Small Magellanic Cloud (SMC; 0.9 dex) (krnault et al.lll988h . The dust- 
to-ga s ratio in Cloud 2 w as estimated to be as low as 0.001 from a submm continuum observa- 
tion (jRuffle et al.ll2007l ). This is about 10 times smaller than in the solar neighborhood, and 
is consistent with the metallicity. Because the Cloud 2 clusters are much closer (D ~ 12kpc) 
than LMC/SMC (D ~ 50kpc), they are quite suitable for the study of t he disk fraction in 



a low-metallicity environmen t. Their youth (most likely age of ~0.5 Myr) (jYasui et al.l 12006 



20081 ; IKobayashi et al.l 120081 ) adds an advantage to the present study because any variation 
in disk fraction should be most clearly seen when the disk fraction value is near its maximum. 
Because the EOG is located at the large heliocentric distance of ~10kpc, high sensitivity 
and high spatial resolution is necessary for the present study. 



Deep JHKs-band images of Cloud 2 clusters were obtained on September 2, 2006 with 
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an exposure time of ~600 sec for e ach band with the 8.2 -m Subaru telescope equipped with a 
wide-field NIR camera, MOIRCS (jlchikawa et al.ll2006l ). The seeing was as good as ~0.35"- 
0.45" through the night. As a result, we obtaine d deeper images than our preliminary studies 
of the Cloud 2 clusters (jYasui et al.lbood . I^OOsh . All the data were reduced with IRAfU with 
standard procedures, and the results are shown in Fig. [TJ JHK aperture photometry has 
been performed using IRAF apphot. The Mauna Kea Observatorie s (MKO) filters were used 
jTokunaga et alJl2002h . GSPC P330-E, which is a MKO standard jLeggett et al.lhoOfih . was 
used for the photometric calibration. Because the fields of the Cloud 2 clusters are very 
crowded, we used the aperture diameters ~0.6" and 0.5" for the Cloud 2-N and 2-S cluster 
members, respectively, to avoid the contamination of adjacent stars. The limiting magnitudes 
(10 a) for the Cloud 2-N cluster are J = 22.3 mag, H = 21.7mag, Kg = 21.2 mag, while 
those for the Cloud 2-S clusters are J = 22.2 mag, H = 21.3 mag, K = 21.0 mag. The 
achieved limiting magnitudes correspond to the mass detection limit of ~0.1M Q , which is 
similar to the mass detection limit of extensive observations of embedded clusters in the 
solar neighborhood with smaller (2-4 m class) telescopes. Therefore, we can compare the 
disk fractions in the low-metallicity environment to those in the solar neighborhood on the 
same basis for the first time. 



3. IDENTIFICATION OF CLUSTER MEMBERS 



Because the Cloud 2 clusters are located behind the Cloud 2 molecular cloud ( Yasui et al. 



20061 ). the cluster members can be clearly identified as extinct red sources behind the CO 
dense cores on the sky (see Fig. [T| see also Fig. 1 in Kobayashi et al. 2008). We identified 
cluster members of the Cloud 2-N and 2-S clusters with the following criteria: (1) distributed 
in regions of the Cloud 2 clusters, which are defined by the yellow box and circle in Fig. (TJ 
and (2) large Ay excess compared to normal field stars (extinction of Ay > 3.5; see Fig. 
2). On the color-magnitude diagram they are located in the region apart from the region 
of normal field stars because the Cloud 2 clusters are located at the peak of the Cloud 2 
molecular cloud core (Kobayashi et al. 2008; Yasui et al. 2008) and reddened by the large 
extinction. We confirmed that stars with Ay > 3.5 mag are concentrated on the cluster 
fields, while stars with Ay < 3.5 mag are widely distributed over the observed field. In Fig. [2] 
the identified cluster members in the cluster regions are shown with red circles while all the 
sources outside the cluster regions (~3.5'x4' region around each cluster) are shown with 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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black dots. As a result, we identified 52 and 59 cluster members for the Cloud 2-N and 2-S 
clusters, respectively, with the 10 a limiting magnitude (red circles in Fig. 2) in both Kg and 
H bands. 

Considering the very large R g of the Cloud 2 clusters (R g ~ 19kpc), contamination 
of background stars is negligible and most of the field objects are foreground stars (some 
are background galaxies). To quantify the contamination by foreground stars, we compared 
the Ay distributions of all the sources in the cluster regions and the field objects outside 
the cluster regions (Fig. [3]). Because the number of field objects decreases significantly at 
Ay ~ 3.5 mag, most cluster members can be distinguished from the field objects as red 
sources with Ay > 3.5 mag. The contamination by the foreground stars are estimated by 
counting the field objects in the tail of the distribution at Ay > 3.5 mag and divide it with 
the total number of the sources in the cluster regions: ~20 % and ~2 % for the Cloud 2-N and 
2-S clusters, respectively. Note that the contamination is very small for the S cluster because 
of the smaller area of the cluster region. Although there must be some cluster members at 
Ay < 3.5 that missed our identification, the number of such sources should not be significant 
in view of the geometry of the molecular clouds, which are in front of the clusters (Yasui et 
al. 2006, 2008). 



4. DISK FRACTION 



We derived disk fractions of the Cloud 2 clusters using a J-Hvs. H- K color-color diagram 
for the identified cluster members (Fig.Hj). Stars with c ircumstellar dust disks are known to 
show a large H — K color excess (ILada fc Adams! Il992l ). On the JHK color-color diagram, 
stars without circumstellar disks are seen as main-sequence stars reddened with extinction, 
while stars with circumstellar disks are seen in "the disk excess region" (to the right of the 
dot-dashed line in (Fig. HJ) because of thermal emission from hot dust disk with a temperature 
of 1000-1500 K. For the Cloud 2-N and 2-S clusters, the disk fractions are 9+4 % (5/52) and 
27±7 % (16/59), respectively, where the ratio of stars with excess to the total number of 
stars is given in parentheses. The uncertainty of the disk fractions is estimated only from the 
Poisson errors because phot ometric errors a re quite small with the median errors of 0.06 mag 
for (H — K) and (J — H) (ILiu et al.ll2003l ). Given the possible systematic uncertainty, we 
cannot say at this stage if the difference between the disk fractions of the Cloud 2-N and 
2-S is significant. Because the above-discussed contamination from the foreground objects 
of Ay > 3.5 mag may decrease the disk fraction, the actual disk fractions could be higher 
than the derived values here. However, the amount of the increase is estimated as +2.2% 
and +0.5% for the Cloud 2-N and -S clusters, respectively, when reducing the total number 
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of the identified cluster members (52 and 59 for the N and S clusters, respectively) by the 
estimated number of contaminated foreground objects (~10 and ~1 for the N and S clusters, 
respectively). Because this maximum possible increase of the disk fraction is well within the 
large Poisson uncertainty (4% and 7% for the Cloud 2-N and 2-S clusters, respectively), 
we conclude that the effect of the background contamination on the disk fraction is small. 
In addition, in the case where we assume that all stars in the cluster regions are cluster 
members, the disk fractions for the Cloud 2-N and 2-S cluster become 6% (5/89) and 21% 
(16/78), respectively. Even in this extreme case, the amounts of the decrease (—3% for the 
Cloud 2-N cluster and —6 % for the 2-S cluster) are again within the Poisson errors. 

To confirm the negligible effect of the photometric errors, we estimated disk fractions of 
two extreme cases: one case is that (J—H) and (H—K) colors of all stars were systematically 
changed within the photometric uncertainties of each source to get the lowest disk fraction 
and the other case to get the highest disk fraction. Assuming a Gaussian distribution of 
the photometric errors, the probability that the color of a star is within the ±1 a offset is 
68 % and the probability that the color of a star is redder than the 1 a offset is about 16 % 
(= (100 — 68)/2). Because the Cloud 2 clusters have about 50 members, the probability 
that all stars have redder colors than the la offset is 0.16 50 , which is essentially zero. If 
we assume that the probability distribution of disk fraction is Gaussian, the estimated disk 
fraction of the extreme case (probability 0.16 50 ) should have a probability less than the 3 
sigma case (0.1%). In the above extreme case (all the sources have +lcr or —la redder 
colors), we found that the disk fraction increases from 9 % to 19 % (N cluster) and from 27% 
to 44% (S cluster). Since such an extreme case only occurs with a probability of 0.16 50 , 
much less than a 3 sigma occurrence, we consider that the 1 a photometric uncertainty must 
be much less than 3 % ((19 - 9)/3 ~ 3) and 6 % ((44 - 27)/3 ~ 6) for the Cloud 2-N and 
2-S clusters, respectively, which is smaller than the Poisson errors. 

In view of the low-met allicity environment of the EOG, there might be some fundamental 
differences between the clusters in the solar neighborhood and in the EOG, such as the 
reddening and the stellar mass distributions, that may require different treatments of the 
disk fraction estimation. In our earlier paper, we confirmed that the reddening vector as 
well as the main-sequence track on the color-color diagram does not change with metallicity 
down to the metallicity of the Cloud 2 (Yasui et al. 2008, section 4). Also, the stellar mass 
distributions of the Cloud clusters are suggested to follow the universal IMF, which is widely 
observed in the solar neighborhood (Yasui et al. 2006, 2008). Therefore, we conclude that the 
use of the same division of with/without disk as for the solar neighborhood is appropriate. 
As for the treatment of the different filter systems, see Yasui et al. (2008) for the detail. 
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JHK DISK FRACTION OF THE NEARBY EMBEDDED CLUSTERS 



We derived the JHK disk fraction of various embedded clusters in the solar neighbor- 
hood using the photometry data in the literature (Table 1 and Fig. [5]). We chose publications 
with the following criteria: i) JHK photometry data of cluster members are available in or- 
der to e stimate disk fraction s in the same way as for our data, ii) the mass detection limit is 
<1 M Q (IHaisch et al.ll2001al ). iii) all data in each embedded cluster are estimated in the same 
filter system. For all clus ters we used the classical dwarf track in the Johnson-Glass sys- 
tem (IBessel fc Brettlll988l ) and a reddening vector in the Arizona system (IRieke fc Lebofsky 
19851 ). JHK disk fractions were estimated in the same way as for the Cloud 2 clusters: 
the border line, which determines whether a star has a disk or doesn't have a disk, passes 
through the point of M6 type stars in the dwarf track on the color-color diagram. 

Recent disk fraction studies often use JHKL data because thermal infrared (A > 3 /mi) 
can trace the colder (thus, outer) disk, and disk fractions from JHKL data are thought 
to be a more robust signature of the circumstellar disk than those only from JHK data 
( IHaisch et al.l l2001aj ) . However, disk fractions from only JHK data are about 0.6 of those 
from JHKL data with slightly more scatter on the disk fraction- age plot, and the estimated 
disk lifetime from JHK data is basically identical to that from JHKL data. Therefore, 
despite a little larger uncertainty, the disk fraction from only JHK data should still be 
effective even without thermal IR data. Since the detection sensitivity in the thermal infrared 
is still not enough for the distant clusters like those in the EOG, we can make use of the 
high detection sensitivity at JHK wavelengths to estimate the disk lifetime for the EOG 
clusters. 



6. DISCUSSION 

6.1. Disk Fraction in Low-metallicity Environment 

We show the disk fractions of the Cloud 2 clusters, as well as of the derived disk fraction 
of seven embedded clusters in the solar neighborhood, in the disk fraction-age diagram 
(Fig. [5j see §EJ). The disk fractions of the Cloud 2-N and 2-S clusters are significantly lower 
than nearby embedded clusters of similar age. As discussed in Section HI the contamination 
from the foreground objects and the possible unidentified cluster members have only a small 
effect on the estimated disk fraction. Although the difference between the filter system for 
our data (MKO) and the data in the literature might cause some systematic offset for the 
disk fraction value, we confirmed that the disk fraction does not depend on the choice of 
the filter system by confirming that the disk fraction of a good reference cluster (NGC 2024 
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cluster) is identical for both CIT and MKO systems (see Appendix A). 



Given the importance of the age of the clusters to the above conclusion, here we discuss 
the uncertainties of the age and the age spread of the clusters. Considering the IMF derived 
from K-band luminosity function (KLF), the age of the Cloud 2 clusters should be less 
than 2Myr: the KLF fi tting shows an unrealistic top-heavy IMF with the age of 2Myr 
( jYasui et al.ll2006l . 120081 1 . Therefore, we conservatively estimate that the age upper-limit of 
the Cloud 2 clusters is l Myr. Moreover, the su ggested supernovae triggered star formation 
of the Cloud 2 clusters (jKobayashi et al.ll2008l ) suggests the age of 0.4 Myr, which strongly 
supports the above estimate from the KLF. Because there is no spectroscopic study of the 
clust er members, the r e is n o observed age spread information as for the nearby clusters as in 
e.g.. |Palla fc Stahlerl (120001 ) . Because the mean age of this cluster is estimated as ~0.5 Myr, 
the age spread should be smaller down to ±0.5 Myr as assumed in the young cluster model 
in Muench et al. (2000). Muench et al. (2002) estimated the age of the Trapezium cluster at 
0.8±0.6Myr. Palla & Stahler (2000) actually obtained smaller age spread for the youngest 
embedded clusters (rho Oph, ONC). Moreover, in the above triggered star formation picture 
(Kobayashi et al. 2008), the age spread of the Cloud 2 clusters should be small because of 
the single triggering mechanism at a certain time. Therefore, we conclude that the most 
likely age of the clusters is 0.5 Myr and no more than 1 Myr. Even if we consider the possible 
age range, the estimated disk fractions are still quite low compared with those in the solar 
neighborhood. 

Although the Cloud 2 clusters appear to be quite close on the sky, they are totally 
independent clusters with a p rojected distance of ~25pc at the heliocentric distance of 
D ~ 12kpc (jYasui et al.l 120081 1 . If they were located at the same heliocentric distance as 
the Orion dark cloud (D ~ 400 pc), the angular separation between the clusters becomes 
3.5°, which is larger than 2.5° between two independent clusters, e.g., NGC 2 024 and NGC 
2071 clusters, listed in the previous study of JHK disk fraction (jLadalll9991 ). Because the 
measured disk fractions for two independent clusters were similarly low, the disk fractions 
of embedded clusters in the low-metallicity environment are inferred to be universally low, 
at least in the EOG. In combination with the fact that there are no very young (<1 Myr) 
embedded clusters of solar metallicity whose disk fractions are this low, our results strongly 
suggest that disk fraction depends strongly on metallicity. 



6.2. Why Lower Disk Fraction with Lower Metallicity? 



The most simple interpretations of the low disk fraction are that th e disk is optically thin 
either because of the collisional agglomeration of dust grains in the disk (jAndrews fc Williams 
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20051 ) or because of the low dust-to-gas ratio due to the low metallicity. However, the for- 
mer is unlikely because the probability of dust collision, which is proportional to the square 
of dust density, should be rather low in the low-metallicity environment. The latter is 
also unlikely because the inner disk should be optically thick even with a metallicity of 
— 1 dex: opacity of the disk with the wavelength of 2/im [y ~ 1.5 x 10 14 [Hz]) at 0.1 AU 
is estimated to be ~10 5 , assuming a mi nimum mass solar nebula (surface gas density of 
~5xl0 4 g/cm 2 at 0.1AU) ( jHayashilll981l) and a standard power-law opacity coefficient k v 
(~15; k v = 0.1 (^/10 12 [Hz]) cm 2 g " 1 ) (iBeckwith et al.lll990l ). and considering the dust-to-gas 



ratio of 0.001 flRuffle et al 



20071 ). which is 1/10 smaller than in the solar metallicity. If we 



assume that the typical disk-to-star mass ratio (0.01; e.g., Andrews & Williams 2005, Natta 
et al. 2006), we would expect the total disk mass of 10~ 3 M Q for the lowest-mass stars in our 
sample (M* ~ 0.1 M©). Because this is roughly one-tenth that of the minimum mass solar 
nebulae (0.01-0.02 M Q ), we still expect a very large optical thickness of ~10 4 at K-band, 
assuming that disk mass measured in the submm correlates directly with surface density in 
the in ner regions. This assumpt ion is used for disks in the solar metallicity environment 



ve-g- 



Andrews fc Williams! 120051 ). Even if we assume the lowest observed disk-to-star mass 
ratio in the nearby star forming regions (~10 -4 ; Andrews & Williams 2005), the expected 
optical thickness is still ~10 2 . Model spectral energy distributions with various disk mass do 
not show a significantly differen t H — K color exc ess (less than 0.05 mag) even if the amount 
of dust disks decreased to 1/10 JWood et alJl2002h . which correspond to the dust-to-gas ratio 
of ~0.001 in the EOG. Therefore, we conclude that the small disk fraction implies that the 
inner region of the disk is cleared out. 

Note that optically thickness does not necessarily mean the existence of observable 
H — K excess for lowest-mass stars wi th low effective temp erature (~3000K) and with very 
low disk mass of 10~ 3 -10 _8 M Q (e.g., lErcolano et al.ll2009l ). Because the lowest mass stars 
of our samples (~O.1M ) are M 6-type stars with an effective temperature (T cS ) of ~3000K 



(see e.g., Fig. 7 bottom figure in lLada &: Ladall2003l ). it may be inevitable that some lowest- 
mass stars in our sample do not have detectable H — K excess simply because of the less 
dust content under low-metallicity, thus contributing to the low disk fraction artificially. To 
investigate this effect, we made a plot of the disk fraction of the Cloud 2 clusters for all 
sources brighter than certain i^-band magnitudes (Fig. [6]). It is clearly seen that the disk 
fractions stay almost flat at the low level and do not decrease even including the M stars 
of the faintest magnitudes (K ~ 20) that may lower the disk fraction as described above. 
Therefore, we conclude that the disk fraction is genuinely low for all sources down to the 
mass detection limit and that clearing out of the inner region of the disk is the reason for it. 
For the above arguments, we note two related st udies. One is the detailed disk fraction study 
of the Trapezium cluster by lLada et al.l (120041 ). who show that even for the latest M-type 
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stars (e.g., M9-type stars, -0.02 M , T cfT ~ 2600 K at 1 Myr) has detectable H — K excess 
and the JHK disk fraction does not change even including the stars of this low-mass. In 
view of the youth of the Cloud 2 clusters (~0.5 Myr) that is similar to the Trapezium cluster, 
we would naively expect similar detectable K-band excess even for the lowest mass stars in 
our sample (~0.1 M , M6-type st a rs) w ith T eff ~ 3000 K even with the low metallicity. Also, 
in the simulation by IWood et al.1 (120021 ). a detectable H — K excess is seen even with very 
small disk mass of 10 -6 M under solar metallicity (which corresponds to 10 -5 M Q with the 
metallicity of Cloud2), assuming a central source with 0.6 M & and T e g = 4000 K. 

Why does the inner region of the disk becomes cleared out effectively with lower metal- 
l icity? There are two well- known mechanisms that can clear the inn er disk: dust sublimation 
( iMillan-Gabet et al.ll2007l ) and the X-wind model ( IShu et al.lll994l ). In the following, we dis- 
cuss if these mechanisms can effectively work with lower metallicity. i) Dust sublimation: 
Since dust sublimation typically occurs around 1500 K, the very inner part of the disk with 
T >1500K becomes dust-free. The dust sublimation radius, _R S , is described approximately 
as (L*/47rT r ^ m cr) 1 / 2 , where is stellar luminosity, T rim (~ 1500 K ) is disk temperature a t 
inner disk radius (R s ), and a is the Stefan-Boltzmann constant (iDullemond et al.l 120071 ). 
However, because the inner disk is expected to be optically thick even with the low metallic- 
ity of —1 dex, the temperature distribution inside the disk for —1 dex should not differ from 
that for the solar metallicity (Odex). Therefore, R s is not expected to be significantly larger 
in the lower-metallicity environment, ii) X-wind model: Due to the interaction between the 
stellar magnetic field with the accretion disk , the gas disk is thought to be truncated at an 
inner radius of R x = ®~£ r ' ' {^JGM^M^) 1 / 7 flshu et al.lfioooh . where M, is the star's mass, 
Md is the disk accretion rate, /i* is magnetic dipole moment, <3>dx is a dimensionless number 
of order unity, and G is the gravitational constant. Among these parameters determining i? x , 
only the disk accretion rate, Md depends on metallicity. The magnetorot ational instability 
(MRI), which is the most likely cause for disk accretion (Hartmann 2009), largely depends 
on the ionization o f the gas disk by Galac tic cosmic rays (jGammie 1996^. X-rays fro m the 
central young star (jlgea fc Glassgold 1999 ). and thermal ionization ( Sano et al. 2000 ). In a 
low- metallicity environment, we can expect that the ionization fraction increases and that 
the dead zone, which is the MRI inactive region of the disk (Gammie 1996), shrinks be- 
cause the recombi nation on g r ain su rfaces is reduced. As a result, Md increases and R x 
decreases. In fact, ISano et al.l (120001 ) suggested that the dead zone shrinks in the case that 
the dust sedimentation of dust grain proceeds, which mimics the case in the low-metallicity. 
On the other hand, low-me tallicity also mea ns fewer molecular ions, which would reduce 



the ionization fraction (e.g.. lSano et al.l 120001 ) . and thus would reduce Md. Which process 



dominates in the low metallicity environment has not been determined conclusively since the 
combination of many physical processes determine the size of the dead zone (e.g, Fromang, 
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Terquem, & Balbus 2002). However, the former process, thus larger Mp, may dominate (T. 
Sano, T. Suzuki, F. Shu, private communication). 

Although further theoretical studies are necessary to have a firm conclusion, we would 
naively expect that R x is not significantly larger in a lower- metallicity environment, and it 
appears that there is no theoretical reason to suspect that lower metallicity would make disks 
less detectable through the H-K excess. The lower dust level would still result in an optically 
thick inner edge and the size of the cleared out inner region should not change significantly 
with the metallicity. Therefore, our observations can simply be interpreted that the inner 
disk (gas and dust) is mostly absent in Cloud 2-N and 2-S. 



6.3. Short Disk Lifetime in a Low-metallicity Environment 

The above arguments lead to the suggestion that the initial disk fraction ( age = in 
Fig. [3P in the low-metallicity environment of our targets is expected to be as high as that for 
the solar metallicity case, assuming that the distribution of initial disk properties is exactly 
the same regardless of metallicity. Although this assumption is not obvious and more detailed 
observations of the disks in the extreme outer Galaxy are needed, it is not an unreasonable 
assumption because the star formation process from gas to stars do not seem to significantly 
change for the metallicity of the Cloud 2 clusters, for example, in view of the universal IMF 
throughout various physical/chemical environments (e.g., Elmegreen et al. 2008, see Omukai 
et al. 2000 for theoretical study). Therefore, we suggest that the significant decrease of the 
disk fraction in Fig. [5] is a result of the much shorter lifetime (~lMyr) of the inner disk 
in the low-metallicity environment com pared to th at in the solar metallicity environment, 
5-6 Myr from JHK disk fraction (§ EJ rLadalll999h . This can be expressed with the rapid 
decrease of the disk fraction as shown with the red thick line in Fig. [51 

We know that the entire circumstellar disk, both the inner and outer disks, dissipates 
almost simultaneously in the solar metallicity environment because th e submillimeter detec- 



tion f raction and the fraction of objects with a NIR excess are identical ([Andrews fc Williams 



20051 ). and also because the observed fraction of transition objects, between the Class II and 



Class III states, is very low (IDuvert et al.ll2000l ). Therefore, we suggest that the entire disk 



dispersal is also rapid in the low-metallicity environment. The distributions of (H—K) excess 
in the Cloud 2 clusters, which are consistent with those of more evolved (~5 Myr) embedded 
clusters with solar metallicity, also support the above idea, though other mechanisms might 
cause the formation of the similar (H—K) distributions. We constructed intrinsic H — K 
color distributions (Fig. [7j) for Cloud 2 clusters and nearby embedded clusters (Table 1). 
The intrinsic (H — K) colors of each star were estimated by dereddening along the reddening 
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vector to the young star locus in the color-color diagram (see Fig. H|). For convenience the 
young star locus was approximated by the extension of the CTTS locus (gray lines in Fig. HJ), 
and only stars that are above the CTTS l ocus were used. F or the Cloud 2 clusters (Fig. [TJ 
top), the CTTS locus in the MKO system (jYasui et al.ll2008l ) was used, while for the nearby 
clusters (Fig. [3, bottom), the original CTTS locus in the CIT system (IMeyer et al.lll997l ) 
was used. To make the figu re clearer, we show the colo r distribution s for only NGC 202 4 
jrlaisch et al.lhoool ). Taurus kenvon fc HartmanrJfl99~5l ). and IC 348 jrlaisch et al.lhoOlbh . 
which have significantly different disk fractions. It is clearly seen that the distribution be- 
comes bluer and sharper with lower disk fractions. The distributions of the Cloud 2 clusters 
resemble that of IC 348, whose disk fraction is as low as the Cloud 2 clusters (~20 ' 



70 . 



6.4. Physical Mechanism for the Rapid Disk Dispersal in a Low-metallicity 

Environment 



As the mechanism of disk dispersal, two possibilities are known : (i) disk accretion 
( IHartmann et al.lll998l ). and (ii) photoevaporation of the surface layer of the gas disk due to 
extreme-ultraviolet (EU V) radiation from the central star and/or the external radiation field 
( iHollenbach et al.ll2000l ). In the low metallicity environment, both possibilities may work 
effectively. First, the disk accretion rate may increase in the low-metallicity as discussed in 
§ 16. 2[ if the increase of the ionization fraction due to less dust is larger than the decrease 
of the ionization fraction due to less molecular ions. This could be the reason for the rapid 
disk dispersal. Secondly, the photoevaporation process may also work effectively to shorten 
the disk lifetime. Because the EUV photoevaporation is solely based on the Stromgren 
condition of the hydrogen gas, metallicity (or dust-to-gas ratio) should not affect the rate of 
the p hotoevaporation mass-loss. However, the recent theoretical work (jGorti fc Hollenbach 
20091 ) shows that far-ultraviolet (FUV) or X-ray radiation is more effective than EUV to 
photoevaporate the disk. Because FUV penetration strongly depends on the amount of dust 
in the disk, they found that the disk lifetime decreases with the dust opacity. Although the 
dust opacity dependence of the disk lifetime is not strong (factor of ~2 changes with the 
dust opacity change of 10), they noted that the results are still in a preliminary stage and 
has to wait further detailed modeling. Therefore, this could be also the reason for the rapid 
disk dispersal in low-metallicit y. In addition, something specific to the EOG environment 



might be the reason (see e.g., lHaywoodll2008l ). but obviously further studies are necessary 



to pin down the physics of this phenomena. 
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6.5. Implication to the Planet-metallicity Correlation 



If the lifetimes of disks at the solar metallicity ([O/H] ~ 0) and at the low-metallicity 
([O/H] ~ —0.7) are approximated as 5Myr and lMyr, respectively, our results suggest that 
the disk lifetime strongly depends on the metallicity (Z = [M/H]) with a ~ 10 z dependence. 
Because the disk lifetime (rdisk) is directly connected to the planet formation probability 
(p p i), the strong metallicity dependence of r disk may create a similarly strong metallicity 
dependence of p p j. Therefore, the disk dispersal could be one of the major driving mecha- 
nisms for the planet-metallicity correlation, which shows a strong sharp (~ 10 2Z ) metallicity 
dependence (Fischer et al. 2005). In the current standard core nucleated accretion model 



of the giant planet formation (e.g., see a review by iLissauer fc Stevenson! 120071 ) . the planet 
formation time scale is determined both by the core growth time and the gas accretion time 
onto the core. Since the first study by Safronov (1969), core growth time in the giant planet 
zone has been thought to be far longer than the observed r disk . But, after the discovery that 
Tdisk is less th an 10 Mvr, it is widely accep ted that giant planet c ores must form via rap id run- 
away growth (IWetherill fc Stewartlll989l ) and oligarchic growth (IKokubo fc Idalll998l ). which 
can grow the cores for the giant planet (5-10 M®) at the outer disk (5AU) in a reasonable 
time, e.g, lOMyr. Because this growth time is comparable to r disk , the variation of Tdi S k due 
to the metallicity around this timescale may cause a strong metallicity dependence of the 
planet formation probability, P p \, thus contributing to the planet-metallicity correlation. 

As far as we know, a deterministic model based on the core-accretion model by Ida & Lin 
(2004) is the only study that quantitatively estimates the effect of Tdisk variation on P p \. They 
show that the planet formation probability decreases with decreasing Tdisk (see their Fig. 2 
(b)). However, the amount of decrease is only a factor of 2-4 with a factor of 10 change in 
the disk lifetime: the dependence is more like P p i ~ r^ k and is not a strong function of Tdisk- 
If this is the case, our suggestion that Tdisk ~ 10 z means t he disk dispersal can contribute 



to the planet-metallicity correlation by a factor of 10 5Z . Ilda fc Lin! (120041 ) show that the 
core-accretion model can qualitatively explain the planet-metallicity relation, but also note 
that the suggested metallicity dependence is not as strong as the observed relation: it is more 
like ~ 10 z dependence judging from their Fig. 2 (b). Therefore, by adding the contribution 
from the disk dispersal (~ 10 a5Z ), most of the observed metallicity dependence of the planet 
formation probability (10 2Z ) can be reasonably explained. However, a quantitative argument 
is difficult at this stage because there are ambiguities and complexities in the planet formation 
models, and thus the planet-formation time scale. Note that other groups are trying to 



reproduce the steep meta llicity dependence by introducing a critical solid mass (IWyatt et al. 
20071 : iMatsuo et al.l 120071 ) for the core accretion model. 



While the planet-metallicity correlation appears to be clear in Fischer's systematic work 
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(IFischer et al.ll2005l ). Udry & Santos suggest that the co rrelation in low-metallicity range 
([Fe/H] < 0) may be relatively weak (jUdry fc Santos! 120071 ) . Santos suggests that the planet 
formation mechanism could be different for the two metallicity ranges: cor e-accretion fo r 
the high-metallicity range and disk instability for the low-metallicity range (jSantod 120081 ) . 
Because our results seem to qualitatively explain a significant part of the 10 2Z dependence 
of the planet-metallicity relation, we assume that the dependence also holds true in the low- 
metallicity range ([Fe/H] < 0), where we made our observations. However, the core-accretion 
model may not be able to explain the existence of planets around the low-metallicity stars 
with such short disk lifetime. 



7. Summary 

We present deep imaging of two low-metallicity star forming regions in the EOG, Digel 
Cloud 2-N and Cloud 2-S. 



1. The fraction of young stars with disks is found to be 9±4 % for Cloud 2-N and 27±7 % 
for Cloud 2-S. Previous works have shown that the age of these clusters is only approx- 
imately 0.5 Myr, and therefore the disk fraction is much lower than expected compared 
to the disk fraction measured for nearby star forming regions. 

2. Various possibilities for the reason for the low disk fraction are investigated. We suggest 
that the simple reduction of the dust-to-gas ratio cannot explain the low disk fraction 
and that the initial disk fraction should be as high as in the solar neighborhood. We 
also suggest that dust sublimation and stellar magnetic activity cannot explain the 
enlargement of the inner disk hole in the low-metallicity environment. Therefore, the 
significant decrease of the disk fraction implies that the disk lifetime is much shorter 
(~lMyr) than in the solar neighborhood (5-6 Myr). We suggest that the short disk 
lifetime is a result of the photoevaporation of the surface layer of the gas disk. 

3. We suggest that the lower disk lifetime with lower metallicity may be a major factor in 
explaining the planet-metallicity correlation. But this is model-dependent and further 
theoretical and observational work is needed. 



Although we derived disk fractions for two very young clusters in this paper, it is neces- 
sary to study more embedded clusters in the EOG to find a more quantitative relationship 
between the disk lifetime and the metallicity. It is also important to study the disks in longer 
wavelengths, e.g., L-band and mid-infrared wavelength, to confirm if disks at larger radii also 
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disappear rapidly, as suggested. These results will be a strong constraint for interpreting the 
planet-metallicity correlation, and ultimately for understan ding how many planet-harboring 
stars we can expect in the Galaxy ( ILineweaver et al.l 12004 1 . 
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A. Comparison of Disk Fractions in Different Filter Systems for the NGC 

2024 Cluster 



NGC 2024 is a good reference target for the comparison of disk fractions in different fil- 
ter systems becaus e it is a very young (0.3 M yr) embedd ed cluster in the solar neighborhood 
( Meyer et al. 1997 ) with a high disk fraction ( Ladal 1999 ). The data were obtained on 1 Jan- 
uary 2007 (based on data collected at the Subaru Telescope and obtained from the SMOKA, 
which is operated by the Astronomy Data Center, National Astronomical Observatory of 
Japan) with a field of view of ~3' x 4' centered on (a 2 oo(b ^2000) — (05 : 41 : 41.4, —01 : 54 : 
59.84). The NGC 2024 data were reduced using the same procedures as the Cloud 2 data. 
For the photometry, we used only stars that are sufficiently isolated with an aperture radius 
of d = 3.4" (~95% energy). For the photometric calibration, we used five stars in the field, 
whose JHK photometric quality in 2MASS is good, after converting 2MASS magnitudes to 
MKO magnitudes. 

Disk fractions can change in different mass ranges and in different spatial regions on the 
sky even when observed in the same filter system. The disk fraction estimated with stars, 
which are in the same region on the sky as the MOIRCS data and in almost the same mass 
range of ~0.7-0.02 M from a previous study in the CIT system (ILevine et al.l 120061 ). was 
55±17% (10/18), which is consistent with the MOIRCS data (65±16%, 17/26) within the 
errors. Moreover, almost all stars in the disk excess region in the CIT color-color diagram 
are located in the same excess region in the MKO color-color diagram. This is true also 
for the stars in the non-disk excess region. Therefore, we conclude that the disk fraction of 
NGC 2024 is correctly estimated in the MKO system. 
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Fig. 1. — NIR pseudocolor images of the Cloud 2 clusters at the Galactic radius of R g ~l9 kpc. 
The color images of Cloud 2-N (top) and 2-S clusters (bottom) are produced by combining the 
J- (1.26 /jm), if- (1.64 /xm), and -fTs-band (2.14 /mi) images obtained at the Subaru telescope 
on September 2008. The field of view of both images is ~3.5'x4'. The yellow box and circle 
mark the locations of the clusters, with closeups shown on the right. 
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Fig. 2.— (J — Kg) vs. Kg color-magnitude diagrams (CMD) of the Cloud 2-N cluster 
(left) and the Cloud 2-S cluster (right). Identified cluster members in the cluster regions 
(yellow box and circle in Fig. [I]) are shown with red circles while all the sources outside 
the cluster regions (~3.5'x4' field for each cluster) are shown with black dots. Only stars 



dwarf tracks ( 


Bessel & Brett 


1988) 


D'Antona & Mazzitelli 


1997. 


1998). 



19881) . while the blue lines show isochrone models (0.5 Myr, 



19981 ). For convenience, the isochrone model was approximated 



by the straight line shown with a solid aqua line. The gray arrows show the reddening 
vectors of Ay = 3.5 mag from the isochrone model, while the dot-dashed aqua lines show 
the isochrone model with the extinction. Stars that are in the Cloud 2 clusters region on 
the sky and are located to the right of the dot-dashed lines on the CMD, are identified as 
the cluster members. The dashed lines show the limiting magnitudes (10 a). 
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Fig. 3. — Ay distributions of all stars in the cluster region (cluster members; black dots and 
line) and all stars outside the cluster region (field objects; gray dots and line). The Ay is 
estimated from Fig [2] based on J — K colors. The distribution of field objects is normalized 
to match with the total area of the cluster regions. Because the number of field objects 
decreases significantly at Ay > 3.5 mag, most cluster members can be distinguished from 
the field objects as red sources with Ay > 3.5 mag. 
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Fig. 4. — (J—H) vs. (H—K) color-color diagrams for the Cloud 2-N cluster (left) and the 
Cloud 2-S cluster (right) in the MKO system. Identified cluster members and field stars are 
shown with red filled circles and black dots, respectively. Typical uncertainty (1 a) of the 
colors are shown at the top-right corner. The solid curve in the lower left portion of each 
diagram is the locus of points corresponding to the unreddened main-sequence stars. The 
dot-dashed line, which intersects the main-sequence curve at the maximum H—Ks values 
(M6 point on the curve) and is parallel to the reddening vector, is the border between stars 
with and without circumstellar disks. The classical T Tauri star (CTTS) loci are shown with 
gray lines. 
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Fig. 5. — Disk fraction as a function of cluster age. JHK disk fractions of the Cloud 2 
clusters are shown with red filled circles, while derived JHK disk fractions of other young 
embedded clusters with solar metallicity are shown with black fi lled circles (see §0). The 
JHKL disk fractions of those other clusters (IHaisch et al.ll2001al ) are shown with gray open 
squares. The black and gray lines show the disk fraction evolution under the solar metallicity, 
while the red line shows the proposed JHK disk fraction evolution in the low-metallicity 
environment. These lines are derived with least-squares fits. 
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Fig. 6. — Disk fraction as a function of the -fT-band manitude. Disk fractions are derived 
for all the sources brighter than the designated magnitude. It is clearly seen that the disk 
fractions stay almost flat at the low level and do not decrease within the uncertainty even 
including the late M stars at the faintest magnitudes (K ~ 20). 
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Fig. 7. — Comparison of intrinsic H — K color distributions. The fraction of stars (/stars) 
per each intrinsic color bin ((H — K)q) are plotted. Top: (H — K)q distributions for the 
Cloud 2 clusters in the MKO system. Bottom: (H — K) n distribut ions for th e various youn g 
clusters in the solar n eighbor hood: NGC 2024 jHaisch et al.lboooh . Cham I jLuhmanlbopj ). 
IC 348 jHaisch et allhoOlbl ). Disk fractions of NGC 2024, Cham I, and IC 348 are 58%, 
41 %, and 21 %. The dashed lines show the borderlines for estimating the disk fraction. Note 
that the border of the (H — K) value for estimating the disk fraction (dashed line) differs 
slightly for different filter system. 
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Table 1: JHK disk fractions of embedded clusters in the solar neighborhood. 



ID 


Cluster 


Age a 


JHK disk fraction 


Mass b 


Filter system 


Reference 01 






(Myr) 


(%) 


(M Q ) 






1 


NGC 2024 


0.3 (1) 


58±7 


0.13 


SQUIID 


2 


2 


Trapezium 


1.5 (3) 


53±3 


-0.03 


SQUIID 


4 


3 


Cham I 


1.8 (3) 


41±4 


-0.03 


CIT 


5 


4 


Taurus 


2.3 (3) 


34±6 


0.3 


SQUIID 


6 


5 


IC 348 


2.3 (3) 


21±4 


0.19 


SQUIID 


7 


6 


NGC 2264 


3.3 (3) 


30±3 


-0.3 


SQUIID 


8 


7 


Tr 37 


4.0 (9) 


17±3 


-0.3 


2MASS 


9 



a Rcfcrences for the ages are shown in the parenthesis. 
b Mass detection limit of the data. 

c The photometric system of the JHK data obtained with the Simultaneous Quad Infrared Imaging De- 



Lada et al. 


19931 


Horner et al. 


1997) 



^References for the JHK photometric data. For Trapezium, onl y the Fred Lawrenc e Whipple Observatory 
(FLWO) data were used to match the mass detection limit (see Muench et al. 20021 ). A Ks filter was used 
instead of K filter for the observation of Tr 37. 
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